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Аннотация: Представлены результаты численного исследования структуры несколько вида 
отрывных течении в канале в зоны их внезапного расширении. Расчеты выполнены на основе 
решения полной системы нестационарных уравнений Навье-Стокса в переменных скорость- 
давление методом разностей на установление. Для численного решения уравнений Навье-Стокса 
использована двух шаговую неявную схему Писмена — Ракфорда. В результате исследованы поля 
скоростей и давления, изучена вихревая структура циркуляционного течении в области за 
уступом и определены протяженность. 


Ключевые слова: уравнения Навье-Стокса, плоский канал с обратным уступом, отрывное 
течение, схему Писмена — Ракфорда, метод контрольного объема. 


1. ВВЕДЕНИЕ 


Исследование отрывных течений имеет важное значение как с фундаментальной, так и 
прикладной точек зрения. Течения жидкости во внезапно расширяющихся каналах встречаются в 
различных технических устройствах и сооружениях. Резкое изменение геометрии стенки канала 
или поверхности обтекаемого тела способно вызвать отрыв потока и существенно изменить его 
кинематическую структуру. Течение в канале с внезапным расширением относится к наиболее 
простому классу отрывных течений, когда точка отрыва потока является фиксированной. 
Теоретический расчет таких течений представляет большие трудности из-за образования сложных 
отрывных и возвратно-циркуляционных течений в области за уступом. Первые расчеты 
стационарных двумерных ламинарных отрывных течений несжимаемой жидкости в каналах были 
изучены Блазиусом еще в 1910 году аналитически в виде рядов [1]. В дальнейшем эта задача 
использовалась многими ученными для изучения механизмов отрывных течений и для 
тестирования разностных схем решения уравнений Навье-Стокса. В силу большой практической 
значимости такие течения изучались теоретически и экспериментально как для ламинарных [2-4], 
так и для турбулентных [5- 7] режимов движения несжимаемой и сжимаемой жидкости. 


Во многих работах в этом направлении рассматриваются течения в каналах с двусторонним 
внезапным расширением [8-14]. Экспериментальные данные для этого случая в канале получены в 
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работах [8, 9], в которых отмечается образование циркуляционной зоны за уступом. Ряд 
исследователей для расчётов течений с внезапным расширением использовали уравнения 
движения в приближении пограничного слоя [14, 15]. В работе [18] было изучена течении 
жидкости в канале с внезапным односторонним расширением который было использована полные 
уравнения Навье-Стокса и взята очень хорошие численные результаты. 


Хорошо известно, что численное решение задач о движении вязкой несжимаемой жидкости на 
основе уравнений Навье-Стокса осложнено не только их нелинейностью, но и отсутствием явного 
уравнения для определения давления. Следует отметить, что многие предыдущие 
экспериментальные исследования проведены для турбулентного движения жидкости, а режим 
ламинарного течения в каналах с внезапным расширением обычно рассматривался при малых 
значениях чисел Рейнольдса но современный прогресс в области компьютерной гидромеханики 
позволяет повысить качество моделирования физических процессов вихреобразования и получить 
устойчивые численные решения при умеренных числах Рейнольдса. В [18] предложен 
эффективный метод численного интегрирования полной системы нестационарных уравнений 
Навье-Стокса в физических переменных скорость-давление для несжимаемой жидкости. 


Цель данной работы состоит в применении указанного метода для решения внутренней задачи 
гидродинамики о отрывные течении жидкости в канале с внезапным односторонним расширением, 
внезапным односторонни дважды расширением, симметричной двухсторонним расширением, 
несимметричной двухсторонним расширением. При этом наряду с изучением структуры течения 
вблизи уступа, преследовалась цель аппробации эффективности данного метода для расчета 
сложных течений, характеризующихся наличием возвратно-вихревых областей течения. 


2. ФИЗИЧЕСКАЯ И МАТЕМАТИЧЕСКАЯ ПОСТАНОВКИ ЗАДАЧИ. 


Рассматривается двумерное течении в канале с различними внезапными расширением в виде 
уступа. Физическая картины анализируемого течении, и конфигурация расчетной области 
представлены на рис. 1. Где а-внезапный односторонний расширения, б-внезапный односторонний 
дважды расширением, в-симметричный двухсторонний расширения, г-несимметричный 
двухсторонний расширения. 
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Рис. 1. Физический схемы течении во внезапно расширяющимся канале 
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Начало введенной декартовой системы координат находится B левом нижнем углу B точке 0. 
Ширина канала в левом входном сечении А, А, имеет размер h, а в правом выходном сечении 


СС, (для а, в на рис.1), СС, (для б, г на рис.1) размер h. Высоты уступов ВВ, (для а на рис.1), 


ВВ, (для б на рис.1) равна соответственноа =й-—й, ВВ,,В,В, (для в на рис.1), BB, D,D, (для г на 
рис.1) равна соответственнор =0.5(Л—,). Левая А, А„ и правая СС, границы расчетной области 
считаются достаточно удаленными от сечения с внезапным расширением BB В. (для а на рис.1), 
В,В,В., ООО, (для б на рис.1), ВВ,,В,В, (для в на рис.1), BB , D,D; (для г на рис.1) чтобы на них 


можно было принять условия, соответствующие невозмущенному течению. Кроме того, 
предполагается, что в канале в области А, А, ВВ, перед уступом выполняются условия полностью 


развитого течения. Тогда в сечении АА, для горизонтальной скорости реализуется 
параболический профиль, форма которого зависит от параметра Е =а/й, Е =Ь/й,. В выходном 
сечении СС, расчетной области вдали от уступа профиль горизонтальной скорости по 
предположению установившийся и описывается параболой Пуазейля. 

Для описания движения жидкости используются нестационарные уравнения Навье-Стокса без 


каких-либо упрощающих предположений. При введении безразмерных величин за масштаб длины 
принимается ширина канала ЁЛ ‚за масштаб скорости принята среднерасходная скорость в канале. 


В безразмерных величинах система нестационарных уравнений Навье-Стокса и уравнения 
неразрывности с постоянными плотностью p = Const и кинематической вязкостью у в декартовых 


координатах имеет следующий вид [23]: 


2 2 
д0 уд уд Фр EREA 
д дх ду дх ôy” ôx 1 
ду ду „ду др ду ду (0 
+U tV += и + | 
д дх ду ду ду” дх 
20 ду _ у 
дх ду 


Здесь U,V — соответственно аксиальная, радиальная составляющие вектора скорости потока, р — 
гидростатическое давление, у - молярная кинематическая вязкость. 

Для завершения постановки задачи необходимо задать начальные и краевые условия на всех 
границах расчетной области. Пусть жидкость движется слева направо для всех каналов. Граничное 
условие для горизонтальной скорости в начальном сечении АА,, а вертикальная скорость У в 
сечении АА, принимается равной нулю. На всех неподвижных твердых стенках выполняются 


очевидные граничные условия прилипания U | г= ди не протекания У ‚= 0, где Г - твердая граница. 


На выходе из расчетной области в сечении СС, для горизонтальной и вертикальной скоростей 


принимаются стандартные условия вытекания Неймана. В начальный момент времени 
горизонтальная скорость в расчетной области имеет соответсвующий параболический профиль, а 
поперечная скорость и давление равны нулю. Основными параметрами задачи являются число 
Рейнольдса и геометрическая высота уступа Е = а/л, Е =Ь/й. Следует подчеркнуть, что давление 
Р в рассматриваемой системе уравнений не является основной переменной ни в одном из этих 
уравнений. При нашем подходе необходимое уравнение для давления выводится из уравнения 
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неразрывности в виде уравнения типа Пуассона. При этом необходимые для его решения значения 
давления в граничных узлах определяются с помощью уравнений движения в комбинации с 
граничными условиями для компонентов скорости [16]. В процессе решения задачи требуется 
определить поля скорости и давления в расчетной области и исследовать влияние числа 
Рейнольдса и геометрического размера уступа Е на структуру течения в канале и протяженность 
циркуляционной области, которая образуется вниз или верх по потоку за уступами. 


3. МЕТОД РЕШЕНИЯ 


Общий принцип решение системы исходных нестационарных уравнений (1) выполняется методом 
конечных разностей на установление. Из-за сложностей согласования полей скорости и давления 
для дискретизации уравнений движения в Х,Ү направлениях и уравнения неразрывности 


использовалась сетка с разнесенной структурой расположения сеточных узлов для зависимых 
переменных. Это означает, что компоненты скоростей и давления определяются в различных 
узлах. Такой подход аналогичен методам SIMPLE [17] и дает определенные преимущества при 
расчете поля давления. Для конечно-разностной аппроксимации исходных уравнений движения 
используются двух шаговую неявную схему Писмена — Ракфорда [19], а для уравнений 
неразрывности используются неявная схема. Схема Писмена — Ракфорда аппроксимации второго 
порядка точности для производных по времени и по пространству т.е. О(лг?,Ах?,лу?). При этом 
диффузионные и конвективных слагаемые аппроксимируются по схеме с центральными 
разностями. Схема Писмена — Ракфорда имеет следующий вид 


На шаге 1 проводится прогонка по направление i 


Ф Ф". + А”. Ф, = Ф, + Лу". Ф ja -Ф,,, — 
1,] tJ LJ 2Ах LJ 2Лу (2) 
п п п Fa кез үс n+l п+1 
= Аг Ф; ы —2Ф,, +Ọ; +At Ф; —2Ф,, + Ф, |, № Ри, T Рл, у | 
Келу? Ве Ах” 2Ах 
а на шаге 2 — прогонка по направление j 
79 — Ф -Ф_, Фф, 
Ф,,-Ф,, + Аг", ИАА дру" БИЕ ЫБ 
” у Е 2Ах у 2 у (3) 
F EN еы %Ћ п+1 п+1 
др iH! 2Ф,, +Ф, Ја д Ф, 2Ф,; +Ф, , L At Pin j Ру 
Re Ay?’ Re Ax’ 2Ах 


Здесь Ф=И приставлен схема только продольного скорости, для поперечного скорости тоже 
использован такой схема. Особенностью дискретизации является то, что конечно-разностная 
аппроксимация центрируется в соответствии с выбранным шаблоном. При этом сеточные индексы 
для зависимых переменных оказываются сдвинутыми [20]. Необходимое уравнение для 
вычисления давления можно получить из уравнения неразрывности. С этой целью будем 
следовать известной процедуре SIMPLE [22]. Отметим, что уравнение Пуасона для давления 
фактически заменяет уравнение неразрывности, и система уравнений оказывается замкнутой. В 
нашем случае эффективным способом решения, рассматриваемого двумерного разностного 
уравнения второго порядка для движении и давлении является его редукция к системам уравнений 
второго порядка с трехдиагональными матрицами, которые решаются методом "прогонки" [20]. 


4. РЕЗУЛЬТАТЫ РАСЧЕТОВ И ИХ ОБСУЖДЕНИЕ 


Данная работа преследовала три основные цели. Первая состояла в апробации численной схемы 
интегрирования двумерных нестационарных уравнений Навье-Стокса в переменных скорость- 
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давление на примере расчета течения в канале с различным внезапным расширением. Вторая - B 
конечно-разностной аппроксимации исходных уравнений использовался двух шаговый неявную 
схема Писмена — Ракфорда. Третьим — в изучении детальной структуры течения и поля скорости 
в канале, на нижней и верхней стенке которого имеется уступы с вертикальным размером 
Е=а/й, Е=Ь/һ. Основные численные расчеты были проведены при значении параметра 


Е=а/й=0.4, E =b/h=0.3 на сетках с шагами Ах = 0.08, Лу = 0.01. Расчеты выполнены для число 
Рейнольдса (Ке= 2000) для каналах внезапным односторонним и односторонни дважды 
расширением. А для симметричной двухсторонним и несимметричной двухсторонним 
расширением расчеты выполнены для число Рейнольдса (Ке=1000). 


На рис. 2 приведены результаты расчетов в виде векторного поля скоростей в расчетной области 
канала внезапным односторонним а) односторонним дважды б) симметричной двухсторонним в) 
несимметричной двухсторонним г) расширением с высотой уступа Е = 0.4 для ї = 10. 
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Рис. 2. Расчетное векторное поле скоростей 


Из приведенных рисунков видно, что для указанных двух чисел Рейнольдса (рис 2), вычисленных 
по выходной ширине канала л, в зоне за уступом характерно образование возвратных течений, а 
протяженность этой зоны и структура циркуляционного течения в ней зависят от числа 
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Рейнольдса и высоты уступа Нетрудно видеть, горизонтальный размер и структура вторичной 
вихрей усложняется и приводит к образованию нескольких локальных вихревых областей, 
которые отчетливо видны (на рис. 2). На рис. 3 приставлен изолиний поле скоростей канала 
внезапным односторонним а) односторонним дважды б) симметричной двухсторонним в) 
несимметричной двухсторонним г) расширением с высотой уступа Е = 0.4 для t=10. 
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E -1--0,9 ш-0,9--0,8 O -0,8--0,7 =O -0,7--0,6 №-0,6--0,5 в-0,5--0,4 Æ -0,4--0,3 
o -0,3--0,2 №-0,2--0,1 Œ -0,1-0 с 0-0,1 E 0,1-0,2 E 0,2-0,3 m 0,3-0,4 
m 0,4-0,5 m 0,5-0,6 m 0,6-0,7 o 0,7-0,8 0,8-0,9 o 0,9-1 E 1-1,1 

Ш 1,1-1,2 E 1,2-1,3 1,3-1,4 Ш 1,4-1,5 1,5-1,6 ш 1,6-1,7 E 1,7-1,8 
ш 1,8-1,9 ш 1,9-2 ш 2-2,1 ш2,1-222 2,2-2,3 ш 2,3-2,4 Ш 2,4-2,5 


Рис. 3. Расчетное изолиний поле скоростей в канале 


Анализ векторного и изолинии поля, и расположение линий равных скоростей, подобных линиям 
тока, наглядно показывает особенности и структуру циркуляционного движения жидкости. Из-за 
прилипания жидкости на нижней и верхней стенке канала движение в циркуляционной зоне 
несимметрично. С ростом числа Рейнольдса центр основного вихря смещается вправо по 
направлению основного потока в канале. Очевидно, что зона рециркуляции в рассмотренных 
случаях подпитывается жидкостью из основного потока благодаря образованию возвратного 
течения в зоне за уступами, как показано на рис. 2-3. 


5. ВЫВОД 


Универсальный дискретный аналог общих нестационарных уравнений Навье-Стокса в 
переменных скорость-давление успешно применен для решения задачи о течении в канале с 
внезапным односторонним расширением, внезапным односторонни дважды расширением, 
симметричной двухсторонним расширением, несимметричной двухсторонним расширением. 
Разработанный численный метод решения системы разностных уравнений движения показал 
эффективность используемой численной двух шаговый неявную схема Писмена — Ракфорда для 
расчета сложных течений с рециркуляциями и высокое качество моделирования тонкой структуры 
течения в вихревых зонах до чисел Рейнольдса Ке < 2000. В результате детального исследования 
полей скорости закономерности формирования вихревой структуры течения в канале с внезапным 
расширением, включая вторичные вихревые образования и нестационарные отрывные режимы 
течения. В заключение следует отметить, что хотя при Ке = 2000 численное решение оставалось 
устойчивым, но относиться к нему следует осторожно и рассматривать его нужно как численный 
эксперимент, полезный для более глубокого понимания свойств решений полных уравнений 
Навье-Стокса. 
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